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ABSTRACT 


A  Dynamic  Simulation  of  Soil-Wheel  Interaction 

Sheridan  L.  Kissklevits 
Advisor 

1.  Robert  Ehrlich 

May  1973 

A  dynaaic  vehicle  Model,  simulating  both  the  vehicle  dynamics  and 
the  wheel-soft  soil  interaction,  is  presented.  Equations  of  notion  are 
presented  for  a  two-dimensional,  six-degree -of-freedoM  vehicle.  Rela¬ 
tionships  for  forces  and  aoasntw  on  the  vehicle  and  its  wheels  due  to 
wheel-soil  interaction  are  developed  aa  functions  of  norsal  and  shear 
stresses  in  the  soil.  Both  the  rigid  wheel  and  flexible  wheel  oases 
are  addressed.  Current  state-of-the-art  relationships  among  shear 
stress,  normal  stress,  wheel  sinkage,  end  wheel  slip  are  used.  Para¬ 
metric  studies  are  used  to  assess  the  qualitative  validity  of  the 
model. 

While  experimental  data  does  not  exist,  comparison  of  the  results 
of  the  simulation  with  known  information  from  the  field  of  off-ths- 
road  Mobility  indicates  the  model  to  be  valid. 
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1 

I.  INTRODUCTION 

The  computer  has  provided  a  valuable  tool  for  use  in  the  planning* 
development,  and  analysis  of  off -the -road  vehicle  systems;  eut-and-try 
techniques  are  giving  way  to  computerised  simulations  as  a  method  of 
an.' lysing  various  systems  in  the  search  for  the  optimum  configuration. 

It  appearu  likely  that  long  procurement  lead  times*  and  their  associ¬ 
ated  skyrocketing  costs,  can  be  reduced  through  the  use  of  simulations. 

A  detailed  dynamic  model  which  includes  vehicle  dynamics  and 
wheel-soft  soil  relationships  i3  necessary  to  realise  the  full  poten¬ 
tial  of  the  computer,  and  to  provide  accurate  modeling  of  prototypes. 
While  they  have  taken  the  necessary  early  steps,  previously  developed 

simulations  are  not  detailed  dynamic  models;  ore  type  depends  upon 

i  2 

equilibrium  equations  *  ,  another  upon  rigid  wheel-rigid  surface  rela- 

3 

tionshipa  ,  while  a  third  models  the  vehicle  components  and  determines 

a 

soil  forces  from  a  data  bank  of  experimental  results  .  This  report  is 
a  significant  deviation  from  most  previous  work  in  that  it  presents  a 
dynamic  vehicle  model  with  particular  emphasis  on  the  wheel-soil  inter¬ 
action.  The  approach  to  the  wheel-soil  problem  considers  the  basic 
shear  and  normal  stresses  acting  at  the  soil-wheel  interface,  and 
derives  all  forest  and  moments  as  functions  of  these  stresses.  Current 
state-of-the-art  relationships  among  shear  stress*  normal  stress, 
sinkage,  -nd  wheel  slip  are  used.  In  this  analysis,  speed  effects  have 
been  neglected,  since  current  literature  indicates  no  velocity  effects 
in  the  range  of  consideration. 
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II.  BACKGROUND 

A.  PREDICTING  WHEEL  PERFORMANCE 

The  validity  of  any  mathematical  simulation  depends  upon  the 
ability  to  formulate  the  relationships  among  the  factors  which  deter¬ 
mine  the  behavior  of  the  prototype.  Of  interest  in  the  sell-wheel 
problem  are  the  relationships  among  shear  and  normal  stresses  within 
the  soil,  forces  on  the  wheel,  and  wheel  slip.  The  formulation  of 
these  relationships  has  proven  to  be  particularly  difficult. 

The  forces  acting  at  the  soil -wheel  interface  are  shown  in 
Figure  1.  Shear  and  norrsal  stresses  each  produce  component  forces  in 
the  horitontal  and  vertical  directions.  The  total  horizontal  and 
vertical  forces  are  the  suss  of  these  components. 

One  of  the  most  widely  accepted  expressions  for  pressure  under  a 
plate  in  soil  is  that  proposed  by  Bekker 


where 

p  -  soil  pressure  on  the  plate 

k  .  k, ,  and  n  are  soil  parameters 

c  f 

B  «  width  of  the  plats 
t  -  depth  of  sinkage  of  the  plate 

Bekker  used  Equation  (1)  to  predict  the  sinkers  of  a  rigid  wheel  in 

soft  soil  from  data  obtained  using  plate  tests. 

Figure  2  shows  a  typical  normal  soil  pressure  distribution  that 
has  been  measured  under  a  rigid  wheel.  Also  shown  is  the  pressure 
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distribution  predicted  by  Bekker'a  equation.  The  significant  differ¬ 
ences  are  evident.  As  can  easily  be  seen,  the  lift  force  predicted  by  using 
Equation  (1)  will  be  larger  than  the  actual  force  realised,  hsnoe  the 

predicted  sinkage  will  be  less  than  the  actual  sinkage. 

6 

Sola  and  Ehrlich  recently  proposed  a  modification  to  Bekker'a 
equation: 

lp‘i’i+^+  (k#+^)‘i  +  ,)n  <2) 

where 

p  -  the  nominal  soil  pressure  under  a  plate 

-  the  initial  soil  bearing  capacity,  independent 
of  plate  shape  (a  parameter  of  the  soil) 

h  «*  the  initial  soil  bearing  capacity,  dependent 
on  plate  shape  (a  parameter  of  the  soil) 

plate  area 

»  hydraulic  radius  of  the  plate  » 

plate  perimeter 

k .  •  the  soil  strength  modulus,  independent  of 
"  plate  shape  (a  parameter  of  the  soil) 

k*  «  the  soil  strength  modulus,  dependent  on 
plate  shape  (a  parameter  of  the  soil) 


1^  -  the  degree  of  soil  compaction,  dependent  on 

previous  loads  on  the  soil  (a  parameter  of  the  soil) 

x  ■  the  sinkage  into  the  soil  at  which  the  nominal 
pressure,  p,  is  measured 

n  "  the  soil  sinkage  exponent  (a  parameter  of  the  soil) 

They  then  used  Equation  (2)  to  predict  the  performance  of  a  rigid  wheel 
in  soft  soil  at  aero  slip  Sela  and  Ehrlich  attempted  to  compensate 


for  the  difference  in  actual  and  predicted  pressure  distributions  by 
neglecting  the  shear  stress  contribution  to  normal  force  on  the  wheel 
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6 

The  performance  pro  dieted  has  shown  good  agreement  with  experimental 
results. 

Q 

The  well-known  Coulomb  equation  gives  the  shear  strength  of  a 

soil: 

r  •  ( c  4*  <rtan0)  (3) 

where 

r  «  shear  strength  (maximum  available 
shear  stress) 

c  *  soil  cohesion 

a  a  normal  pressure  on  the  soil 

jj  *  angle  of  internal  soil  friction 

o 

Using  Equation  (3).  Janosi  and  Hananoto  proposed  the  following  equa¬ 
tion  to  describe  a  soil  shear  stress-strain  curve: 

t  •  (c  +  vtan0)(l  a  e“*^x)  (4) 

where 

a  *  the  base  for  Naperian  logarithms 

j  »  the  soil  def  oraation 

K  *  the  deformation  modulus  of  the  soil 
shear  stress-strain  curve 

As  applied  to  a  loaded  soil-shear  plate,  this  equation  states  that  to 
displace  the  plate  a  distance  j,  the  shear  stress  r  Bust  be  overcome. 

Note  that  this  is  less  than  the  maximum  shear  strength  given  by  Equa¬ 
tion  (3),  and  that  for  very  large  deformations  ( j  -a-  « ) ,  Equation  (4) 
approaches  Coulomb's  equation. 

In  extending  this  equation  to  predict  shear  stress  distribution 
10 

under  a  wheel,  Janosi  assumed  that  the  wheel  contact  patch  was  a 
series  of  small  plates  that  anted  independently,  and  obtained  an 


expression  for  j  froa  the  cycloidal  movement  of  the  wheal  periphery. 

Severe!  other  theories  regarding  shear  stress  have  been  advanced, 
bat  they  fall  into  two  categories: 

(1)  the  solutions  are  incomplete  because  the  number  of 
unknowns  exceeds  the  number  of  equations  which  can  be  written. 

(2)  if  complete,  the  solution  is  limited  to  a  very  special 
soil  condition. 

Accordingly,  J  ariosi*  s  original  approach  is  still  the  one  used  by  many 
researchers. 

It  is  believed  that  sinkage  and  normal  stress  are  functions  of 
wheel  slip.  However,  these  relationships  have  not  yet  been  discovered. 
B.  VEHICLE  SIMULATION 

In  an  extensive  report  on  vehicle-terrain  simulations,  Schuring 
and  Belsdorf  1  have  presented  a  model  of  a  vehicle  traversing  soft, 
smooth  soil.  The  simulation  assumed  steady-state  movement,  and  hence 
was  postulated  in  terms  of  equations  of  equilibrium. 

The  wheel  performance  equations  were  derived  from  an  assumed 
triangular  distribution  of  normal  stress,  with  the  maximum  stress 
determined  by  Bekker's  equation}  from  the  Janosi-Hanaaoto  equations  for 
shear  stress  distribution  (simplified  by  some  assumptions  involving  the 
relative  magnitudes  of  sines  and  cosines);  and  froa  an  assumed  empi¬ 
rical  relationship  between  leading  and  trailing  portions  of  the  contact 
sons.  See  Figures  3  and  h. 

Execution  of  the  simulation  requires  simultaneous  solution  of 
various  algebraic  equations,  and  provides  the  following  output: 
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•  the  "thrust"  force  at  the  contact  rone  of  a  driven  wheel 

•  the  drive  torque  at  the  axle  of  a  driven  wheel 

•  the  slip  of  the  driven  wheel 

•  th9  forward  velocity  of  the  vehicle 

The  model,  however,  has  not  been  verified  by  experiment. 

In  the  report  of  their  study  of  egress  performance  of  vehicles, 

2 

Muddappa  and  Baker  have  presented  a  model  of  a  vehicle  traversing 
soft,  smooth  soil.  The  simulation  determined  inertial  forces  acting  on 
the  vehicle  as  a  result  of  its  velocity,  and  assumed  the  vehicle  and 
forces  acting  on  it  to  be  in  dynamic  equilibrium  (steady-state)  at 
finite  increments  of  time.  The  simulation  was  postulated  in  terms  of  a 
highly  non-linear  set  of  differential  equations. 

The  wheel  performance  equations  were  derived  from  an  assumed 
linear  contact  patch;  from  the  wheel  sinJcage  given  by  Equation  (1);  and 
from  the  "corapactir  resistance"  of  the  soil,  computed  by  integrating 
Equation  (1),  See  Figure  5. 

Execution  of  the  simulation  consists  of  an  iterative  procedure, 

and  provides  the  position  of  the  vehicle  as  part  of  its  output.  The 

model  also  has  not  been  verified  by  experiment. 

4 

Schoeh  and  Shah  ,  in  their  report  on  the  development  of  computer- 
modeling  techniques,  have  presented  an  example  model  of  a  tracked 
vehicle  traversing  a  soil  terrain.  This  simulation  used  extensive 
models  of  vehicle  components  and  their  interaction  to  calculate  engine 

torque  and  speed  of  wheel  sprockets. 

The  track  performance  was  determined  from  a  table  of  test  values 
for  a  eimiliar  vehicle.  Using  sprocket  speeds  as  input,  and 

- yt  -a  iWir  -  - 


considering  factors  such  as  slope,  rolling  resistance,  and  soil  coali¬ 
tions,  the  model  determines  sprocket  forces  as  its  output.  Vehicle 

speed  follows  as  a  function  of  sprocket  forces. 

3 

Sloss,  Ehrlich,  and  Worden  have  presented  a  model  of  a  vehicle 
e grossing  rrcra  water  onto  a  hard  ramp.  The  simulation  assumed  the 
vehicle  to  be  in  non-steady-state  motion,  and  was  postulated  in  terms 
of  its  differential  equations  of  notion.  Since  the  simulation  assumed 
a  hard  bank,  the  wheel  performance  equations  derived  used  normal  forces 
and  the  coefficient  of  friction  to  determine  wheel-bank  forces.  Execu¬ 
tion  of  the  simulation  involved  solving  the  differential  equations  of 
notion  by  a  numerical  method,  with  outputs  of  vehicle  location,  speed, 
acceleration,  and  pitch  attitude.  The  computer  model  was  verified  by 
scale  node!  tests  and  has  shown  dose  agreement  with  experimental 


results 
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III.  THE  MATHEMATICAL  MODEL 


A.  THE  VEHICLE  AND  ITS  REFERENCE  FRAMES 

The  simulatio:.  discussed  in  this  report  uses  the  US  Army  M151A1 
£  Ton  Truck,  or  "Jeep*,  as  its  prototype  vehicle.  The  vehicle  is 
assumed  to  be  two-dimensional  (vertical  plane)  with  6  degrees  of  free- 
don.  The  hull  of  the  vehicle  has  freedom  of  surge,  heave,  and  pitch*. 
Each  wheel  is  free  to  heave.  In  addition,  the  front  and  rear  wheels 
are  free  to  rotate,  but  they  are  not  independent  because  of  the  action 
of  the  transfer  case. 

There  are  three  coordinate  systems  used  in  the  simulation;  they 
are  shown  in  Figure  6.  The  X*Z'  system  is  the  inertial  reference 
frame,  with  the  X*  axis  lying  along  the  horizontal.  The  XZ  system  is 
affixed  to  the  vehicle,  with  its  origin  at  the  vehicle  center  of 
gravity.  When  the  vehicle  is  at  rest  on  a  horizontal  profile,  the  X 
exis  is  horizontal.  The  xz  system  is  affixed  to  the  hub  of  each  wheel, 
and  carries  a  subscript  which  denotes  the  specific  wheel  which  it  de¬ 
scribes  (1  -  Front,  2  -  Rear).  The  xz  axes  are  parallel  to  the  respec¬ 
tive  XZ  axes.  Note  that  since  the  wheel  is  only  free  to  heave,  the  x 
axis  will  not  be  used  in  this  model. 

The  vehicle  suspension  is  simulated  by  a  spring  and  damper  system. 
Since  there  is  generally  very  little  movement  of  wheels  in  the  horizon¬ 
tal  direction,  the  suspension  is  assuasd  to  be  solidly  affixed  to  the 


Surge  is  forward  motion;  heave  i3  up-and-down  motion;  and  pitch 
is  rotation  about  an  fj  a  which  goe3  through  the  center  of  gravity  of 
the  vehicle,  and  is  perpendicular  to  the  plane  of  the  vehicle. 
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vehicle  with  no  freedom  of  movement  in  the  X  direction. 

B.  EQUATIONS  OF  MOTION 

To  represent  properly  the  acceleration,  velocity t  anu  location  of 
a  vehicle  undergoing  movement,  we  must  rely  upon  time-dependent  equa¬ 
tions  of  motion.  These  equations  for  a  rigid  body  have  been  derived 
’1 

many  times  before  ,  so  their  derivation  is  omitted  here.  Instead, 
they  are  presented  in  the  generally  accepted  form,  with  comments  as  to 
how  they  are  applied  to  the  model. 

The  equations  used  in  this  simulation  are  the  minimum  required  for 
program  debugging  and  initial  model  validation.  In  reality,  wheels 
accelerating  with  respect  to  the  body  will  produce  inertial  forces  on 
the  body.  These  forces  have  been  neglected  since  this  is  a  first-cut 
model  designed  to  perform  on  uniformly  sloping  soil  with  uniform 
strength  (hence,  nearly  equal  wheel,  sink  ages) .  Under  these  circum¬ 
stances,  the  inertial  forces  mentioned  are  small.  They  must  be  in¬ 
cluded  later,  if,  for  example,  the  model  is  used  to  study  obstacle 
climbing. 

1.  Hull  Equations 

Six  degrees  of  freedom  of  a  solid  body  give  rise  to  six 
equations  of  motion.  Of  these  six,  the  equations  for  surge,  heave,  and 
pitch  are: 

Surge;  Fx  -  M(U  +  QW  -  RV)  (5) 

Heave:  F.  -  H(W  +  FV  -  QU)  (6) 

Pitch:  Gy  -  IyQ  -  P(I#R  -  I^P)  ♦  Ed,?  -  ImR)  (7) 

where 

D,7,W  -  velocities  in  ths  X,  Y,  and  Z  directions 
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P,Q,R  «  rates  of  rotations  about  the  X,T,  and 
Z  axes 

M  -  the  mass  of  the  body 

FX,F„  •  the  net  forces  on  the  body  in  the  X 
and  Z  directions 

I  -  the  product  of  inertia  of  the  body 
about  the  X  and  Z  axes 

G  -  the  net  torque  on  the  body  about  the 
7  X  axis 

«  aonents  of  inertia  of  the  body  about 
7  the  X,Y,  and  Z  axes 

Under  the  two-dime naional  assumption,  the  terms  V,P,  and  E  are  sero,  so 


the  equations  reduce  to  the  following! 


Fx  -  H(U  +  QW) 
F_  -  M(W  -  QU) 

z 

G  -IQ 

y  y 


(8) 

(9) 

(10) 


2.  Wheel  Equations 

The  two  degrees  of  freedom  at  the  wheels  give  rise  to  two 
equations  of  notion:  for  heave  and  for  pitch  (Equations  (6)  and  (7)). 
Again,  the  terns  V,P,  and  R  are  sero,  so  these  equations  reduce  to 

Fsi  “  Ki(^i  '  QiUi>  1  -  1,2  (il) 

G  -IQ  (12) 

yv  yww 

where  the  subscripts  denote  front  and  rear  wheels  according  to  the 

notation  of  Figure  6,  and 

G  -  the  combined  net  torque  on  both  wheels 

yv 

»  the  combined  wheel-drive  train  rotary 
7  moment  of  inertia 

Q  -  Q.  -  Q2  (since  the  wheels  are  not 
independent) 
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C.  SOIL  FORCES  ON  A  RIGID  WHEEL 

In  order  to  get  expressions  for  the  forces  end  moments  in  the 
equations  of  notion,  I  now  turn  to  an  aralysis  of  the  wheel-soil  inter¬ 
action. 

All  forces  exerted  on  the  wheel  by  the  soil  are  assuned  to  arise 
froa  two  basic  stresses:  normal  stress  and  shear  stress.  The  effects 
of  these  stresses  are  to  resist  sinkage  and  forward  notion  of  the  wheel 
due  to  compaction  on  the  soil,  to  resist  rotation  of  the  wheel  due  to 
shearing  of  the  soil,  and  to  provide  traction.  Each  stress  will  be 
considered  individually,  and  examined  for  its  affect  ou  wheel  motion. 
The  rigid  wheel  case  will  be  treated  first,  since  it  is  applicable  over 
a  considerable  range  of  sinkages  for  certain  soils,  and  it  is  the 
simplest  case  to  discuss. 

1.  Normal  Stress  Forces  and  Hoaents  On  a  Rigid  Wheel 

The  normal  stress  is  directed  radially  inward  on  the  wheel. 

A  compressive  stress  is  considered  positive,  no  negative  stresses  are 
allowed,  and  the  stress  is  assumed  to  be  uniform  across  the  width  of 
the  wheel.  The  stress  at  a  given  point  is  the  average  stress  of  an 
elemental  force  acting  on  an  elemental  area.  Points  on  the  wheel  sur¬ 
face  are  located  by  an  angle  y  measured  froa  the  horizontal.  A 
counterclockwise  angle  is  positive.  The  angle  to  the  point  where  the 
wheel  first  touches  the  soil  ("^)  is  called  the  "entrance  angle*,  and 
the  angle  to  the  point  where  the  wheel  leaves  the  coil  C^)  ie  called 

the  *-2xit  angle."  The  area  of  soil -wheel  contact  between  these  two 
angles  is  called  the  "contact  patch."  Referring  to  Figure  ?,  the 


>y//X\W///// 


FIGURE  ?  NORMAL  STRESS  UNDER  A  RIGID  WHEEL 


normal  stress  is  expressed  as 


dK 

a  » 

b  r  &y 

where  b  is  the  width  of  the  wheel. 


(13) 


The  normal  stress  produces  forces  which  can  be  resolved  into  com¬ 
ponent s  in  the  vehicle  coordinate  system.  This  is  shown  in  Figure  7. 
The  eoaponent  of  the  force  on  the  wheel  in  the  X  direction  due  to 
noraal  stress  is  given  as 

dFXff  »  -  dH  sin  a  (14) 

Siwllinrly,  the  component  of  the  foroe  on  the  wheel  in  the  Z  direction 


due  to  normal  stress  is 

dF„„  -  -  dN  cosa  (15) 

oU 

From  Figure  7 

a  •  -j-  +  7  -  8  •  -j-  +  (7-6)  (16) 

Fro*  trigonometry 

-  cos(7  -  0)  (17) 

■  -  sin (7  -  0)  (18) 

(19) 

sod  (15)  yields 
dFXff  "-bra  cos (7  -  0)d7  (20) 

dfXff  »  b  r  tj  sin(v  -  9)d7  (21) 

To  get  the  total  force  on  the  wheel  in  each  direction.  Equations  (20) 
and  (21)  are  integrated  across  the  entire  contact  patch,  which  is 


sin 


C08 


[t  +  (7-8>] 
’>] 


f+<T- 


Fro*  Equation  (13) 

dN  *  a  b  r  d7 

Substituting  these  into  Equations  (14) 


Siailiarly,  the  component  of  the  force  on  the  wheel  in  the  Z  direction 
doe  to  shear  stress  is 

d Ftr  -  -  dT  sina  (30) 

From  Equation  (28) 

dT  »  r  b  r  d7  (31) 

Making  this  substitution,  and  the  same  substitution  fora  as  before, 
Equations  (29)  and  (3^)  become 

d?XT  »  -  b  r  r  sin^  -  0)d7  (32) 

dFSr  «  -  b  r  r  cos(7  -  0)d7  (33) 

Again,  integration  yields  the  total  force  on  the  wheel  in  each  direc¬ 
tion: 


XT 


r  sin(7  -  9)d7 

w 

r2 

r7’ 

T  cos(7  -  0)d7 

(35) 

Since  the  shear  stress  acts  tangentially,  it  tends  to  retard  tht 
rotation  of  the  wheel.  This  can  be  measured  as  a  moment  at  the  wheel 
center.  Designating  this  moment  Mr,  and  referring  again  to  Figure  9, 

dMT  -  r  dT  (36) 

Making  the  previous  substitution  for  dT,  and  integrating  across  the 
contact  patch  gives 


Mr  »  b  r 


If*1' 

J:*' 


(37) 


As  in  the  case  of  normal  stress,  the  forces  due  to  shear  stress 
produce  a  moment  about  the  vehicle  center  of  gravity.  Again,  only  the 
component  force  in  the  X  direction  is  transmitted  directly  to  the  body. 


The  development  of  the  expression  for  is  precisely  the  sene  as  the 
development  of  M^,  with  ths  obvious  substitution  of  d?XT  for  dF^, 
Hence, 

2  fy' 

Mxt  »  (Z^  +  21)Fxt  +  b  r  t  ain(7  -  0)coa(7  -  0)d7  (38) 

Equations  (34-),  (35).  (37).  and  (38)  give  the  forces  and  nonsnts 
which  result  froa  shear  stresses  on  the  wheel.  Figure  9  shows  that  Fyr 
resists  sinkage  of  the  wheel,  hence  supports  the  vehicle,  while  FJr 
tends  to  move  the  wheel  in  the  forward  direction;  hence  provides 
traction. 

D.  NORMAL  AND  SHEAR  STRESS  FUNCTIONS 

The  task  now  is  to  express  a  and  r  as  functions  of  7,  Numerous 
12  13 

studies  *  have  shown  that  normal  stress  is  generally  distributed  in 


or,  in  stapler  fora. 


Y  (Degrees) 


FIGURE  10  ASSUMED  NGRMAL  STRESS  DISTRIBUTION 


FIGURE  11  LOCAL  GROUND  SLOPE  AND  SINKAGE  FOR  A  RIGID  FRONT  WHEEL 
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idler* 


o  -  A77  +  B7  +  C 


(7,  -  72  )(72  -  7,  ) 

B  -  -  (7,  +  7a)A 
C  «  7,72A 

The  maximum  normal  stress,  <7m„  ,  occurs  at  the  center  of  the  con¬ 
tact  patch.  Bekker’s  equation  for  pressure.  Equation  (1),  is  used  to 
determine  <rm,„  .  This  requires  an  expression  for  the  depth  of  sinkage 
at  the  center  of  the  contact  patch.  Figure  IS  shows  a  rigid  wheel  and 
a  general  soil  profile.  Since  the  coil  is  generally  not  exactly  hori¬ 
zontal,  or  uniformly  sloping,  the  angle  %  is  called  the  "local  ground 
slope."  It  is  the  angle  between  the  horizontal  and  the  chord  connec¬ 
ting  the  points  of  intersection  of  the  wheel  with  the  undisturbed  soil 
profile. at  the  front  of  the  wheel,  and  the  original  soil  profile  at  the 
rear  of  the  wheel.  All  sinkages  are  measured  perpendicular  to  this 
chord.  The  sinkage  is  given  as 


The  maximum  pressure  is  ncy  given  as 


sB  -  r  sxn(0G  -  7H)  -  r  ein(8^  -  7,) 
-  r[»l”<8G  -  \>  -  sin(90  -  >,)] 


T -N 


One©  the  front  wheels  pass  a  point  on  the  soil  profile,  the  soil 
is  compacted,  and  the  rear  wheels  do  not  see  the  "original"  soil  pro¬ 
file.  The  normal  stress  distribution  will  now  be  modified  to  account 


. . 1 
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for  this  prscompaction. 


Figaro  12  shows  typical  experimental  pressor*  vs  sinkage  curves 
14 

presented  by  Bekker  .  These  corves  show  how  the  pressure  acts  when 


the  load  on  a  plate  is  relieved,  then  re-applied.  The  curves  show  that 


when  the  load  is  re-applied,  the  pressure  almost  immediately  returns  to 


its  value  at  the  tins  the  load  was  relieved  from  the  plate.  Soil  engi¬ 


neers  have  used  this  behavior  pattern  for  some  tine  to  aid  in  the  con¬ 


struction  of  Void  Ratio  vs  log  Pressure  curves  used  in  studying  con¬ 


solidation  settlements  ,  It  is  also  the  basis  of  the  tens  in  the 


Sela-Ehrlich  equation  for  pressure  (Equation  (2)).  This  pressure  re¬ 


bound  is  accounted  for  as  shown  in  Figure  13.  The  normal  stress  dis¬ 


tribution  curve  is  as  stated  to  be  between  the  points  located  by  the 


theoretical  entrance  angle  with  the  original  soil  profile,  called  71o 


and  the  actual  exit  angle  with  the  current,  or  existing,  soil  profile. 


The  distribution  curve  is  not  extended  to  the  theoretical  exit  angle 


with  the  original  soil  profile,  72o,  since  the  pressure  vs  sinkage 


curves  show  that  the  pressure  alaost  immediately  drops  to  xero  when  the 


load  is  relieved.  With  these  modifications ,  the  faotors  in  Equation 


(41)  become,  for  the  rear  wheel, 


<>,*  -y2)(y2  -  y,J 


b  -  -  ( y,0  +  72)ji 


c  -  71o72A 


The  sinkage  at  the  point  of  maximum  stress,  which  is  now  at  the  center 


of  the  "modified*  contact  patch,  is 


FIGURE  13  NORMAL  STRESS  DISTRIBUTION  UNDER  A  REAR  WHEEL 
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For  simplicity.  Equations  (45)  -  (48)  will  be  used  for  both  front  and 
rear  wheels,  since  for  the  front  wheel,  7to  is  just  equal  to  7, . 

As  noted  in  Section  IIA,  there  exist  several  theories  on  the 
distribution  of  shear  stress  under  a  wheel  undergoing  slip.  .  Janoai’s 
approach  is  employed  sonoe  it  is  generally  used.  Schuring  and 
Belsdorf  1  presented  an  expression  for  j.  As  modified  for  the  conven¬ 
tions  and  variables  used  in  this  report,  the  expression  is 


J  -  r 


|V,  -  7  +  (1  -  s)^COS(7  -  6)  -  COs(7,  -  e)^  j 


where  a  is  the  wheel  slip  given  by 


in  which 

r  »  the  radius  of  the  wheel 

aw  *  the  angular  velocity  of  the  wheel 

i  «  the  velocity  of  the  vehicle 

Substitution  of  Equation  (49)  into  Equation  (4)  gives  the  expression 


for  shear  stress: 


r  »  (c  +  crtan^)  I  1  -  exp 


(4[-  - 


-  a)^coa(7  -  8)  -  cos (7,  -  8) 
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S.  INTEGRATION  OF  THE  SOIL  FORCE  AND  MOMENT  EQUATIONS  FOR  THE  RIGID 
WHEEL 

1,  Integrations  Involving  Normal  Stress 

Equation  (41)  makes  it  possible  to  integrate  directly  Equa¬ 
tions  (22),  (23),  and  (2?)  for  F_  ,  F_  ,  and  .  With  the  substitu- 

a<T  mo  X” 

tion  for  <x  ,  these  equations  beooas 

f  Yl 

Fx<r  -  -  b  r /(A^  +  BY  +  C)cos(Y  -  d)  dY  (52 


F  »br 

tor 


-Yi 

•/W 

J'ir, 


+  BY  +  C)sin(Y  -  6)  dY 


Mx*  -  <ZRi  +  Zi*x<r  + 


/-Vx 

b  r2 1  {k  f 
•/v 


+  BY  +  C)sin(Y  -  9)cos(Y  -  0)  dY 


Those  are  involved,  but  straightforward  problems  in  integration  by 
parts.  Ths  results  are: 


(Y2  -  Yic)cos(Ys  -  0)  + 


FXff  -  b  r  |  -  (A 'l  +  BVi  ♦  C)sin(Yx  -  0)  +  A  j(Y2  -  Yic)cos(y2  -  0)  + 
(>a  V/x  ♦  Yxc)cos(Yi  -  0)  +  2(sin(Yi  -  0)  -  sin(Y2  -  0))j|  (55) 

a  •  b  r  |-  (AY^  +  BYx  +  C)cos(Yi  -  0)  +  A^Y^  -  Y  )sin(Ys  -  0)  + 
^Yj  —  Y10  "  -  ®)  ♦  2 (cOB\  Y^  -  9)  -  C0#(X>  (56) 


F*c  -  b  r 


)sin(Y  -  0)  + 

3  3 


(2Y  Y 
i  —  xo 
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H**  -  (Zjii  +  Zi)FXff  4.  b  r2  |  (a/  +  BY  +  Cttsin2(Y.  -  6)  + 
4l2(Y  •  Y  )(Y  -  Y  )  +  (Y  -  Y  )sin(2(Y  -  0))  + 

8J^  10X13  10  3  '  3  7 


(2Y  -  Y  +  Y  )ain(2(Y  -  ©))  ♦  coa(2(Y  - 

110  3  '  1  '  '  1 


008 


(2(V  -  8))]| 


(57) 


When  evaluated  fa?  a  given  wheel  position.  Equations  (55)  and  (57) 
are  inserted  into  the  hull  equations  of  notion,  while  Equation  (56)  is 
inserted  into  the  wheel  equations  of  notion  as  shown  below. 

2.  Integrations  Involving  Shear  Stress 

Equation  (51)  makes  it  extremely  difficult  to  integrate 
directly  Equations  04),  (35),  (3?),  and  (38)  for  Fxr,  ?tff  HT,  and 
A  numerical  integration  simplifies  the  problem  considerably,  and  is 
easily  inpleiaented  on  a  computer.  An  accurate  method  is  the  Newton- 
Cotes  5-Point  Method.  As  an  illustration,  the  integral  of  the  shear 
stress  equation  was  calculated  by  two  methods  for  various  values  of  Y^ 

and  v  .  One  method  was  the  Newton-Cotes  Method,  the  other  was  the 

3 

Trapeaoid  Method,  using  intervals  of  five  degrees.  The  Trapezoid 
Method  is  known  to  be  an  accurate  numerical  method  for  computing  inte¬ 
grals,  but  it  requires  many  mors  computations  than  the  Newton-Cotea 
Method,  Tbs  results  are  shown  in  Tabic  1.  In  all  cases,  the  value 
computed  by  the  Newton-Cotes  Method  was  within  O.kSi  of  the  value 
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TABLE  1 

COMPARISON  0?  NUMERICAL  INTEGRATION  METHODS 


SLIP 

*1  ^ 

AREA  BY  NEWTON  -COTES 

AREA  BY  TRAPEZOID 

DEVIATION 

OSC 

°o 

H 

1 

°o 

1 

4.461263 

4.481281 

0.45j6 

-30°  -110° 

2.516806 

2.529037 

0.4036 

c % 

-40°  -110° 

1.156235 

1.161734 

0.4756 

2036 

-20°  -110° 

5.170 466 

5.187433 

0.3# 

20f> 

-30°  -1100 

3.327870 

3.33783? 

0.3036 

zoi 

-40°  -110° 

1.890615 

1.894482 

0.2036 

4036 

-20°  -110° 

5.543893 

5.562296 

0.33* 

4 Oi 

-30°  -iio° 

3.759944 

3.771866 

0.3236 

40$ 

-40°  -110° 

2.294688 

2.300244 

0.2456 

6o£ 

-20°  -110° 

5.757924 

5.777705 

0.3456 

60f^ 

-30°  -110° 

4.007252 

4.021091 

0.3456 

605C 

-40°  -110° 

2.530277 

2.527608 

0.2936 

8036 

-20°  -110° 

5.890605 

5.910926 

0.346 

803b 

-30°  -110° 

4.158725 

4.17372? 

G.3# 

803& 

Jf0°  -110° 

2.675515 

2.684094 

0.3236 

loogt 

-20°  -110° 

5.978328 

5.998789 

0.3436 

10056 

0  O 

-30  -110 

4.257007 

4.272829 

0.3756 

10056 

-40°  -110° 

2.769586 

2.779135 

0.3436 
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computed  by  the  Trapexoid  Method.  When  evaluated,  these  force  and 
moment  terms  are  inserted  into  the  equations  of  notion  in  a  Manner 
siadliar  to  that  for  normal  stress. 

F.  SOIL  FORCES  ON  A  FLEXIBLE  WHEEL 
1.  Tire  Deflection 

As  long  as  the  pressure  exerted  on  a  tire  by  a  soil  stays  be¬ 
low  sons  maximum  value,  the  tire  acts  like  a  rigid  wheel.  When  the 
soil  pressure  reaches  this  maximum  value,  the  tire  deflects.  The  maxi- 
mas  allowable  value  of  soil  pressure  under  a  wheel  is  assumed  to  be  ex¬ 
pressible  as 

P'-P1  +  Pe  (58) 

where 

i 

P  “  the  maximum  soil  pressure 
?j«*  the  inflation  pressure  within  the  tire 
PQ*  the  carcass  pressure  of  the  tire 
This  assumption,  and  the  assumption  of  a  synnetricai  normal  stress 
distribution  with  maximum  value  at  the  centerline  of  the  modified  con¬ 
tact  patch  (Section  III.D)  leads  to  the  following  analysis: 

a.  The  first  deflection  of  the  tire  will  be  at  the 
centerline  of  the  modified  contact  patch  if  the  maximum  soil  pressure 
reaches  p' 

b.  Since  the  soil  pressure  cannot  exceed  p'  the  normal 
pressure  distribution  may  be  based  upon  a  theoretical  sinkago  which 
would  obtain  if  the  wheel  remained  rigid.  This  theoretical  sink  ago  is 
called  and  the  corresponding  theoretical  maximum  normal  stress  is 
called  ^rnixTH* 
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c.  As  ffm,xTH  increases,  more  of  the  tire  deflects.  The 
“tone  of  tire  deflection*  is  defined  by  that  portion  of  the  theoretical 
normal  pressure  distribution  which  has  ordinates  greater  than  K  Note 
that  the  sons  of  tire  deflection  is  symmetrical  about  the  centerline  of 
the  modified  contact  patch,  and  that  the  normal  pressure  in  the  cone 
is  p! 

Additionally,  the  tire  deflection  is  assumed  to  be  linear.  A  typically 
deflected  tire  and  its  associated  normal  stress  distribution  as  assumed, 
in  this  report  is  shown  in  Figure  14. 

Justification  for  thie  assumed  stress  distribution  has  been  pre¬ 
sented  by  Ellis  ^ ,  who  shows  that  the  pressure  under  a  flexible  tire 
on  a  hard  surface  is  nearly  uniform  across  the  contact  patch,  and  drops 
off  sharply  at  the  leading  and  trailing  edges.  Since  the  soil  pressure 
ie  at  the  maximum  allowable  under  the  deflected  portion  of  the  tire, 
the  coil  will  not  deflect  any  further,  hence  it  acts  -3  a  rigid  surface. 
The  remainder  oX  the  tire  still  in  contact  with  the  soil  will  be  hare 
according  to  the  rigid  wheel  theory  described  above. 

The  problem  of  integrating  Equations  (22),  (23),  and  (27)  is 
cite  red  only  by  the  need  t-o  know  two  additional  angles  ( Y  and  Y  )  for 
limits  of  integration.  For  the  flexible  wheel  case.  Equation  (22),  for 


r-1689 


FIGURE  14  NORMAL  STRESS  DISTRIBUTION 
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Consider  the  equation  for  norsaal  stress  distribution, 

c  «  K?  +  ?/  +C  (60) 

It  is  necessary  to  solve  this  for  the  angles  when  o  »  Note  that  in 
the  expression  for  A,  ffm,x  *  jmixTH  .  Substituting  for  a  in  Equation 
(60)  givss 

Ay2  +  By  +  C  -  p'  (61) 

or 

AY*  +  BY  +  (C  -  ?)  -  0  (62) 

Applying  the  quadratic  formula  to  solve  for  y  gives 


components  in  the  vehicle  coordinate  system.  This  angle  is  the  angle 
between  the  centerline  of  the  deflection  tone,  and  the  Z  axis.  From 
Figure  15*  the  angle  is 


«'«  +  Y  )  -  9  (6?) 

10  3 

Note  also  that  the  expression  for  dN  (Equation  (19))  as  no  longer  valid 
since  the  incremental  area  is  no  longer  equal  to  b  r  dv.  However,  it 
is  not  necessary  to  integrate  actually  any  functions  in  the  zone  of 
deflection.  Because  of  the  symmetry  and  constant  pressure  in  the  zone 
of  deflection,  the  resultant  of  all  stresses  in  the  zone  acts  at  the 
centerline  of  the  modified  contact  patch.  To  find  the  resultant, 
consider  the  contact  patch  in  the  zone  of  deflection.  Since  it  is  the 
chord  of  a  circle,  the  length  of  the  patch  is  given  by 


L  .  2  r  sin(«vi(  .  V  „)) 

(68) 

!  f. 

i  ,3 

"si 

The  resultant  is  than 

'4 

F  -2brP  aixifti  V,  -  V  )) 

(69) 

'i; 

z 

X 

The  components  of  the  resultant  force  in  the  X  and  Z  directions 

are: 

5r 

I 

•  Fcz*x“  -  2  b  r  P  ■*»(*<  -  Vsp»  sina1 

> 

¥ 

-  -  2  b  r  P'  sin(KVi?  -  Y^})  cos(i(Yi{.  +  Yg)  -  e) 

(70) 

i 

Fota^  -  **  2  b  r  Psin(*(Ylp  -  Ya  ))  cooa* 

' 

»  2  b  r  ?'  sin (K  v  -  Y  ))  8in(}(  y  +  Y  )  -  0) 

X  IP  JP  '  X  10  3  ' 

(71) 

Equations  (22)  and  (23)  can  now  be  modified  for  the  flexible  wheel: 


*iSaaaBfl& 


F!C  -  -  b  r 


P  fa  Coa(Y  -  9)  dY  -  b  T  fa  coa(Y  -  8)  dY 

“  V  */>. 


X  'Y 

2  ip 


->  i. 


b  r  p'  sJ.tt(±<Yip  -  V  P  ))cos(i(Yi0  +  Ys  )  -  0) 
Y  Y 

/*2P  /*  1 

«  b  r /<rsin(Y  -  8)  dY  +  b  rla  sin(Y  -  8)  dY 

2  \p 

+  2  b  r  Pain(KYif)  -  Y, })ain(i(Yio  ♦  Y^ )  -  e) 


(72) 


(73) 


Evan  though  the  stress  in  the  sons  of  deflection  no  longer  acta 
radially,  it  still  does  not  produce  a  moment  on  the  wheel  because  it  is 
symmetrical  about  the  centerline  of  the  contact  patch.  However,  the 
moment  about  the  vehicle  CG  due  to  the  components  of  forces  due  to  nor¬ 
mal  stress  is  affected,  because  the  moment  arm  is  no  longer  measured 
from  the  edge  of  a  circle,  as  it  was  for  the  rigid  wbc-el. 

Since  the  resultant  F..,,  acts  at  the  center  of  the  contact  Bone, 
the  moment  arw  from  its  X  component  to  the  x  axis  of  the  wheel  is 

d  »  r  cos(KYip  -  Ygp ))cos( at’)  (74) 

The  ’foment  about  the  vehicle  CG  due  to  the  X  component  of  Fctg  is 

«’•  +h  +  d)  <«> 


M*«  -  2  b  r  P*  sini(v  -  Y  )cos(£<y  ♦  Y  )  -  0)[z„.  + 

IP  2P  '  10  2  ^  ftl 

Z,  -  r  cos±{y  -  y  )ain(i(Y  +  y  )  -  0)j 

*  IP  2P  '  10  2  'J 

Finally,  the  expression  for  follows  from  Equation  (24) 

Y  Y 


(76) 


*Xa 


/*2  P 

•fh i  +  zi 


+  r  coaa)  dF  +  M 

X0 


'+  f(ZRi  +  Zi  +  r  cosa)  dFXff  (77) 

♦A# 


zp 
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In  order  to  be  able  to  write  this  in  a  somewhat  condensed  font,  let 

V  -  i<V  -  Y  )  ( 

16  p  ip  sp 


0  -  $<Y  +  V  )  -  0 

10  a 


Then  Y  y 

Hz.  -  b  r  <Z[ti  +  Zl>  W  co8(v  -  0)  dy  -  J^r  cos(y  -  0)  dyj  - 

2  IP 

2  b  r  P*sin(v^  )cos(i8)|zRi  +  -  r  eosiy^  )sin(d)j  + 

r  y  yj 

b  r2|^rsin(V  -  0)cos(Y  -  )  dY  + ^ain(Y  -  0)cos(Y  -  8)dY|  (80) 
*-  a  ip  -* 

Equations  (72),  (73),  and  (80)  give  the  forces  and  moments  which 

result  from  normal  stresses  on  a  flexible  wheel. 

3.  Shear  Stress  Forces  and  Moments  On  a  Flexible  Wheel 

Figure  16  shows  the  shear  stresses  acting  on  a  deflected  tint. 
As  in  the  case  of  normal  stress,  at  any  point  along  the  deflection  sons, 
the  same  angle  is  used  to  resolve  the  forces  due  to  shear  stress  into 
components  in  the  vehicle  coordinate  system.  The  approach  to  finding 
the  resultant  of  the  shear  stresses  acting  in  the  deflection  kobo  is 
the  same  one  used  for  normal  stress.  It  does,  however,  require  the 
simplifying  assumption  that  while  She  shear  stress  distribution  is  not 
syranetrical  about  the  centerline  of  the  contact  patch,  the  eccentricity 
is  small  enough  to  neglect. 

Figure  17  shows  a  typical  shear  stress  distribution  under  a  de¬ 
flected  tire  as  calculated  from  Equation  (51),  using  a  ■  p'in  the  sone 
of  deflection.  Note  that  it  is  not  syranetrical  about  the  centerline  of 
t-he  contact  patch.  The  centroid  of  the  area  bounded  by  Y  and  Y 

1  P  QP 


H  -  b  r 

X<r 


2  b  r  P  sin(y  )cos(|8) 

‘/tf 


6  f 


If  ; 
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(and  hence,  the  point  of  action  of  the  resultant),  does  not  lie  on  the 
centerline  of  the  tone  of  deflection.  The  eccentricity  in  terse  of 


is  approxiaately  1  .  This  eccentricity  affects  two  values:  the  aoaent 


ar*  of  F„„  about  the  wheel  center,  and  the  moaent  am  froa  F_  to 
Car  Car* 


the  x  axis  of  the  wheel.  Both  of  these  values  are  affected  because  of 


the  tom  r  cos{y,  ).  For  values  of  Y,  which  can  reasonably  be  ex- 

J$P  YjV 


pected  (Y  <  20  ),  the  difference  between  r  cos(Y  )  and  r  cos(Y  ♦  l  ) 

A?  ht  ;>p 

is  negligible. 

The  resultant  of  the  forces  due  to  shear  stress  in  the  zone  of 
deflection  can  now  be  written  as 


F _ «2br  sin(Y 

oer 


i 


IP 

rdY 


(81) 


2P 


The  coaponenta  of  the  resultant  force  in  the  X  and  Z  directions  are: 


Y 


PC*rX  ”  2  b  X*  Sin(V>  C0S“jj  ^ 


2P 


-  -  2  b  r  sin(Y„  )  sin(/S)/r  dY 


IP 


(82) 


-  -  2  b  r  sin(Y,,  )  si nah  dY 
cars  ^p  / 

2  P 


y 

SP 

Yip 


,»Yip 

«  -  2  b  r  sin(Y,ip)  c os05)/r  dY 


(83) 


Equations  (34)  and  (35)  can  now  be  modified  for  the  flexible  wheel: 

/*Ys  p  /*Yi 

l(  Y  -  8)  dY  -  b 

■4:p 


/*Y2  P  /-Yl 

Fxt  •  -  b  r  It  sin(  Y  -  8)  dY  -  b  r  fr  sin(  Y  *  8)  dY 


ip 


-  2  b  r  ain(  )  si n(3)j^  dY 


(84) 
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T  »  -  b  r 


riz  i *  /*Vx 

rfr  co8(Y  -  0)  dY  -  b  r/r  cos(y  -  8)  d\ 


/-Yip 

-  2  b  r  ain(Y^)  cos(/3)/r  dY 


(85) 


The  expression  for  the  roment  on  the  wheel  of  force®  doe  to  sheer 
stress,  Equation  (37) »  is  easily  modified  for  the  flexible  wheel  by 
using  the  resultant  F  and  its  moment  arm,  r  coe(Yi/  ). 

ClT  /JP 


*Ya  p 

Mt  «  b  r^/r  dY  +  2  b  r2 


if 


/*YlP  ,  /*Yl 

sin(Y^  )  cos(v^  )/r  dY  +  b  r d/  (86) 

Q  P  T1 » 


The  SKHoent  about  the  mehicle  OG  due  to  the  X  component  of  F(^r  is 

M"- W*R1  + vd>  W> 


H"«-2br  ain(Y^e )  sinf/S)^^  +  +■ 

Yxf 


r  Cos(Y^? )  sin(i3)J  £  dY 
J*Aip 


(88) 


Then  the  expression  for  MXr  follows  from  Equation  (24): 

■Yfep 


+  Zi  +  r  cosa)  d?xr  (89) 


r  ts  p  /•  •* 

Mxr  -/ (ZRi  *  2i  *  r  «**“>  ^xr  *  H"  +  Ki 
Hfe  *\ip 

"  b  r  ^i  +  V  [•i  •in(v  -  9)  dY  -  £«  -  »)  *rj- 
2  b  r  sinty.^ )  sia{0){z^  ♦  -  r  coa(Y,_,s )  ain(3)jjCdY 


'fep 


b  rZ  jl  sizZ(y  -  8)  dY'  +^8in2(Y  -  9)  dY' 


(90) 


Equations  (84),  (Spit  (86),  and  (90)  give  the  forces  and  moments 
which  result  £roa  shear  stresses  on  a  flexible  wheel. 
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G.  INTEGRATION  0?  THE  SOIL  FORCE  AND  MOMENT  EQUATIONS  FUR  THE  FLEXIBUS 
WHEEL 

1.  Integrations  Involving  Noraal  S trails 

As  in  the  rigid  vhael  case,  the  expressions  Tor  ?xc,  F  , 
and  Hxa  are  analytically  integraOle .  The  results  of  the  integrations 
frost  the  rigid  wheel  case  are  directly  applicable  to  the  flexible  wheel 
case,  since  the  only  difference  is  the  lixits  of  integration.  The  re¬ 
sults  are: 


Fx<r  -  b  r 


!(Ay  2  +  By 

ip  ip 


+  C)ein(Y  -  0)  -  (AY  2+  By  +  C)sin(Y  -  6) 
ip  ii  i 

r 


F  *  b  r 

s«r 


-  (AY  *  +  BY  +  C)sin(Y  -  0)  +  A  (Y  -  Y  )cos(Y  -  0) 

s*  ap  st  [_  a  is  a 

+  (Y  -  2y  ♦  Y  )cos(y  -  0)  +  2(sin(Y  -  0)  -  ein(v  -  0)) 

♦  ( Y  -  Y  )coa(Y  -  0)  ♦  (Y  -  2Y  +  Y  )oos(y  -  d) 

IP  1C  IP  IP  1  10  1 

♦  2(sin(Yi  -  0)  -  sin(Yi?  -  0))j  -  2  ?'  sin(Y^)  cos(/3)|  (9- 

<(AY  s  tBV 

1  ip  ip 

-  (AY  2  +  BY  +  C)coa(Y  -  0)  +  A  < 

3P  2P  SP  L 

+  (2Y  -  y  -  Y  )ain(Y  -  0)  +  2(cos(y  -  0)  -  coa(Y  -  e)) 

a»  a  _  10  a  p  '  2  p  a 


+  C)oos(Yip  -  0)  -  (AY^+  BY  *  C)sos(Yi  -  0) 


(Y  -  Y  }cos( Y  -  0) 
io  a  a 


♦  (Y  -  Y  )sin(Y  -  0)  +  (2Y  -  Y  -  Y  )sin(  Y  -  0) 

10  IP  IP  1  10  IP  1 

+  2(coa(Yi  -  6)  -  coa(v  p  .  9))J  *  z  P*  sin(Y^?)  ain(0)j 


(92) 
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"  *ZRi  +  +  2  b  )  cos(y^)  sin(0)  cos(0) 

+  b  rZ|(AY  2  +  By  ♦  C)isia2(Y  -  9) 

|  2?  Sf  2P 

+  (Ay  3«*  By  +  C)i«in2(Y  -  9)  -  (Ay  3  +  By  +  C)isin2(Y  -  9) 

i  i  i  ip  ip  '  ip 

-  »r 

+*ft  2(Y  -  Y  )(Y  -  Y  )  +  (Y  -  Y  )ain(2(Y  -  9)) 

O  10  2f  2P  2  10  2  \  2  '' 

■» 

♦  (2\,  +vs).in(2(vf  -e))  4-  COS  (2(Ysf  .  0))-  cos(2(Yg  -9)) 

+  2(Yio  -  Yi)(Yi  -  Yif)  +  (Y  o  -  Yi().in(2{Yif  -  0>) 

+  «Y  *  Yc  +  )sin(2(Yi  -  0))  +  cos(2(Yi  -  9)) 

-  co8(2<Yi?  -  9))jj  (93) 

Whan  evaluated  for  a  given  wheel  position.  Equations  (91)  and  (93) 
are  Inserted  into  the  hull  equations  of  motion,  while  Equation  (92)  is 
inserted  into  the  wheel  equations  of  notion  as  shown  below. 

2.  Integrations  Involving  Shew-  Stress 

As  in  the  rigid  wheel  case,  the  expressions  for  FXt,  Far,  Mr, 
and  MXr  are  not  easily  integrable.  Again,  a  numerical  solution  is  the 
best  approach  in  this  case.  When  evaluated,  these  force  and  moment 
terms  are  inserted  into  the  equations  of  motion  in  a  manner  siadliar  to 
that  for  normal  stress. 

H.  EQUATIONS  OF  MOTION  IN  THE  SIMULATION 

With  the  soil  forces  and  moments  iust  derived,  Equations  (8),  (9), 
(10),  (11),  and  (12)  can  now  be  used  to  specify  the  particular  equa¬ 
tions  of  motion  for  the  system  used  in  the  simulation.  Figure  18  shows 
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the  vehicle  with  the  forces  and  aonenta  acting  on  it.  Froa  Figure  18, 
and  Equations  (8),  (9),  and  (10),  the  hull  equations  are 


jfas  +  ««  -Ifn  *  ?x2  -  W'«io9]  <« 

Z’CO  -  — J-jw'cort  -  D,  -  D2]  (95) 

Iy9  -  DjXj  .  +  Mr,  +  KrJ  -  \t  -  (96) 

where 

and  D2  ***  suspension  forces 
*X  “  *X<r  ♦  ?xr 

-  the  engine  torque  on  a  wheel 
the  subscripts  denote  front  and  rear  wheels 
Figure  18  also  shows  a  uhoel  with  the  forces  and  assents  acting  on 
it.  Froa  Figure  18,  and  Equations  (11)  and  (12),  the  wheel  equations 


where 


h  *  Vco  "iqj*!  °°’e  +  D1  -  ?.l] 
*2  *  Vco  "4[*2  C°S0  *  ®2  ‘  ?«z] 


V«  '  tl  *  Hv2  '  Hr,  -  Mr, 


e  »  the  angular  velocity  of  the  tfceels 
(the  saw  for  all  wheelo) 


F  •  F  4-  F 
*  tff  gtr 


Note  that  the  signs  of  and  K_  are  reversed  froa  the  hull  equations, 
since  for  the  wheel,  clockwise  rotation  is  positive. 


IV.  THE  SIMULATION 


A.  GENERAL 

The  nodal  in  this  roport  was  developed  as  part  of  a  simulation  in 
a  larger  project  headed  by  Dr.  K.  Peter  Jurkat,  Davidson  Laboratory, 
Stevens  Institute  of  Technology.  The  simulation  used  as  a  subroutine 
a  program  which  the  author  wrote  to  determine  the  forces  and  moments  on 
the  vehicle  due  to  the  wheel-soil  interaction. 

B.  CONCEPT 

Unbalanced  forces  are  calculated  as  a  result  of  vehicle  velocity, 
and  wheel  position  with  respect  to  the  vehicle  and  the  ground.  These 
unbalanced  forces  are  used  to  move  the  vehicle  to  a  new  position  and 
velocity  via  the  equations  of  motion. 

At  the  beginning  of  the  simulation,  initial  position  and  velocity 
are  assigned  to  the  vehicle.  This  position  determines  the  amount  of 
wheel-soil  interferance  in  terms  of  the  exit  and  entry  angle  of  the 
wheel.  Using  these  arid  the  soil  parameters,  Equations  (34),  (35),  (37)* 
(38),  (55),  (56),  and  (5?)  calculate  the  forces  and  moments  derived  in 
Section  HI  for  the  rigid  wheel  case.  If  the  flexible  wheel  case  ob¬ 
tains,  Equations  (84),  (85),  (86),  (90),  (91),  (92),  and  (93)  are  used. 
These  forces  and  moments  are  then  inserted  into  the  equations  of  motion, 
which  are  numerically  integrated  to  get  a  new  position,  velocity,  and 
acceleration.  The  process  is  then  repeated. 

C.  DESCRIPTION 

The  simulation  program  consists  of  five  major  subroutines,  which 
are  shown  in  block  diagram  in  Figure  19.  The  simulation  is  quite 


5 


SOIL  FORCES 

Computation  of  Soil 
Forces  and  Moments 


FIGURE  19  PROGRAM  SUBROUTINES 
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flexible  with  respect  to  types  of  vehicles,  soil  parameters,  end  ground 
profiles  which  may  be  used.  A  description  of  the  subroutines,  with 
comments  on  their  capabilities,  follows, 

IIi?UT 

(1)  Reads  vehicle  description;  the  vehiclo  sprung  weight, 
unsprung  weight,  moment  of  inertia,  location  of  CG,  etc.,  must  be 
specified. 

(2)  Reads  soil  profile  and  parameters.  The  program  stores 
the  soil  profile  as  coordinate  locations,  so  any  profile  can  be  used. 

The  soil  parameters  can  be  varied  along  the  profile  to  simulate 
changing  soil  conditions  which  the  prototype  might  encounter. 

(3)  Reads  initial  conditions  of  the  problem.  The  vehicle 
can  be  started  from  rest,  or  with  an  initial  velocity. 

(4)  Positions  vehicle  with  respect  to  the  soil  profile, 

(5)  Reads  integration  control  variables. 

INTEGRATION  CONTPjOL 

(1)  Calls  a  subroutine  which  integrates  the  equations  of 
motion  by  the  Hamming  Predictor-Corrector  Method. 

(2)  Tests  for  completion  of  simulation. 

UPDATE 

(1)  Calculates  wheel  position  with  respect  to  soil  profile 
in  terms  of  y  ,  Y  •  and  V  . 

X  10  3 

(2)  Calculates  local  ground  slope  and  wheel  slip. 

(3)  Calls  soil  force  subroutine. 

(4)  Calculates  engine  torque;  checks  engine  torque  vs  RFH 
map  and  shifts  to  proper  transmission  gear. 


v* .-.-i  ^ 


Since  no  experimental  data  exists  which  relates  vehicle  perfor¬ 
mance  to  soil  parameters,  the  nodal  was  "validated"  on  a  qualitative 
basis.  It  is  generally  expected  that  vehicle  performance  will  increase 
as  soil  strength  increases.  To  determine  if  the  model  did  produce  re¬ 
sults  with  this  trend,  the  following  approach  was  selected: 

•  Part  i 

1.  Select  c,  kQ,  n,  K,  0.  and  a  uniformly  sloping  soil 

profile. 

2.  Vary  to  change  soil  strength, 

3.  Use  the  ability  of  the  vehicle  to  climb  the  slope  as  a 
measure  of  the  effect  of  soil  strength. 

•  Part  2 

For  a  soil  strength  which  allows  the  vehicle  to  just  climb 
the  slope,  change  the  per  cent  of  slope  and  examine  its  effect  on  ve¬ 
hicle  performance. 

For  the  first  portion  of  the  validation,  the  soil  selected  was  a 
cohesionless  one  with  the  following  parameters: 

o  -  0 
kc  -0 
n  -  1 

K  -  0.2 

tan$  -  0.28 

A  1 5%  slope  was  selected.  Values  of  chosen  were  30»  *K>,  50. 

The  higher  values  of  relate  to  stronger  soil,  as  seen  from  Equations 


(44)  and  (51).  The  initial  velocity  selected  was  60  inches  per  second, 
or  3.4-  miles  per  hoor. 

Table  2  shows  the  vehicle  velocity  at  0.1  second  intervals  for 
each  case.  Figure  20  is  a  plot  of  velocity  vs  tine  for  each  case.  The 
first  1.5  seconds  is  considered  the  transient  portion  of  the  simla- 
tion;  iv,  arises  from  the  assignment  of  initial  conditions  which  are 
not  steady-state.  For  a  of  30,  the  vehicle  could  not  continue  to 
climb  the  slope,  and  the  velocity  gradually  decreased.  At  a  kj  of  4$, 
the  vehicle  was  able  to  slightly  accelerate  up  the  slope.  At  a  of 
50,  the  vehicle  was  able  to  accelerate  even  acre.  Thus,  as  soil 
strength  increased,  vehicle  performance  increased  as  expected. 

For  the  second  portion  of  the  validation,  the  value  of  40  for  k^ 
was  selected.  The  simulation  was  executed  for  values  of  slope  of  1?.5£ 
and  20f>.  Table  3  shows  the  vehicle  velocity  at  0.1  second  intervals 
for  each  ease.  The  plot  in  Figure  21  shows  that  as  slope  is  increased, 
vehicle  performance  decreases  as  expected. 
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TABI2  2  \ 


VEHICIE  VELOCITY  AT  SELECTED  TIKES  FOR  VARIOUS 

TIKE  (sec) 

VELOCITY  (in/sec) 

-  30 

■••40 

0 

VA 

1 

0,00 

60.00 

60.00 

60.00 

0,10 

62.14 

60.66 

60.84 

0.20 

60.26 

60.46 

61.18 

0.30 

58.58 

60.89 

62.15 

0.40 

58.60 

61.48 

63.27 

0.50 

58.84 

62.0? 

64.3? 

0.60 

58.90 

62.49 

65.35 

0.70 

58.78 

62.80 

66.09 

0.80 

58.16 

63.07 

66.81 

0.90 

57.47 

63.43 

67.64 

1.00 

56.56 

63.90 

68.46 

1.10 

56.17 

64.31 

69.30 

1.20 

56.01 

64.74 

70.19 

1.30 

55.86 

65.24 

71.40 

1.40 

55.86 

65.93 

71.86 

1.50 

56.55 

65.94 

72.56 

1.60 

56.28 

66.18 

73.27 

1.70 

54.55 

66.38 

73.98 

1.80 

54.03 

66.66 

74.73 

1.90 

53.68 

66.97 

75.48 

2.00 

53.33 

67.2? 

76.21 

2.10 

53.05 

67.59 

76.96 

2.20 

52.84 

67.85 

77.62 

2.30 

52.52 

68.12 

78.34 

2.40 

52.22 

63.38 

78.95 

2.50 

51.91 

68.68 

79.64 

2.60 

5 1.65 

68.96 

80.28 

2.70 

51.32 

69.26 

81.08 

2.80 

51.1C 

69.58 

81.84 

2,90 

50.78 

69.91 

82.52 

3.00 

50.46 

70.21 

83.21 

ELOCITY  VS  TIME  FOR  VARIOUS  k 
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TABLE  3 

VEHICLE  VELOCITY  AT  SELECTED  TIMES  PO? I  VARIOUS  SLOPES 


155f  SLOPS 


l?j£  SLOPE 


2$  SLOPE 


0.00 

60.00 

60.00 

60.00 

e.10 

60.66 

60.53 

59.5? 

0.20 

60.46 

59.09 

56.97 

0.30 

60.89 

58.7? 

55.90 

0.40 

61.48 

58.61 

55.15 

0.50 

62.0? 

58.50 

54.31 

0.60 

62.49 

58.23 

53.38 

0.70 

62.80 

57.89 

52.28 

0.80 

63.0? 

57.48 

51.04 

0.90 

63>3 

57.08 

49.86 

1.00 

63.90 

56.75 

48.78 

1.10 

64.31 

56.47 

47.72 

1.20 

64.74 

56.22 

46.69 

1.30 

65.24 

55.52 

45.65 

1.40 

65.93 

55.18 

44.5? 

1.50 

65.94 

55.65 

43.36 

1.60 

66.18 

54.84 

43.11 

1.70 

66.38 

54.22 

42.76 

1.80 

66.66 

53.80 

40.72 

1.90 

66.97 

53.30 

39.22 

2.00 

67.27 

52.80 

37.95 

2.10 

67.59 

52.61 

36.70 

2.20 

67.85 

51.69 

35.48 

2.30 

68.12 

50.93 

34.20 

2.40 

68.38 

50.56 

33.05 

2.50 

68.68 

50.19 

31.43 

2.60 

68.96 

49.84 

29.1? 

2.70 

69.26 

49.36 

27.10 

2.80 

69.58 

49.14 

25.61 

2.90 

69.91 

48.12 

24.16 

3.00 

70.21 

47.50 

22.56 
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Yl.  CONSLDSIOKS  ASP  RE 


moss 


A.  CONCLUSION 

On  *  qualitative  baaie,  the  model  of  vehicle  dynamics  and  vhael- 
soft  soil  interaction  presented  is  this  study  is  a  valid  simulation  of 
actual  vehicle  performance. 

B.  RECOMMENDATIONS 

1.  That  experimental  data  he  gathered  to  assess  the  validity  of 
the  model  on  a  quantitative  basis. 

2.  That  the  equations  of  motion  presented  in  this  stncty  be  modi¬ 
fied  to  include  the  relative  accelerations  of  the  wheels  discussed  in 
Section  HI.B. 

3.  That  be  soil-wheel  equations  presented  be  modified  to  use  the 


Sela-Ehrlich  improved  equations  for  normal  stress. 


-Kjt.  ^iSwgvSt ^^gzg^-%S%&V5^pi&i3Zi 
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